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Distance-redshift data can impose strong constraints on dark energy models even when the equa- 
tion of state is oscillatory. Despite the double integral dependence of the distance on the equation 
of state, precision measurement of the distance-redshift relation for z = — 2 is more incisive than 
the linear growth factor, CMB last scattering surface distance, and the age of the universe in distin- 
guishing oscillatory behavior from an average behavior. While oscillating models might help solve 
the coincidence problem (since acceleration occurs periodically), next generation observations will 
strongly constrain such possibilities. 



I. INTRODUCTION 

The subtle variations in the cosmic expansion hold fun- 
damental information on the nature of dark energy and 
the cosmological model. While these variations are not so 
extreme as to cause nonmonotonicity in the scale factor 
vs. time relation a(t), and even less so in its integral the 
distance-redshift relation, we might consider nonmono- 
tonicity in the dark energy equation of state w{a). Hav- 
ing more than one period of strongly negative equation 
of state, and hence acceleration, could ameliorate the co- 
incidence problem of why acceleration is happening now 
out of all the expansion history of the universe. Here 
we examine the implications of periodicity (e.g. a sin- 
gle Fourier mode of nonmonotonicity) in the equation of 
state w{a), and what constraints can be placed on it. 

First, we note that many periodic or nonmonotonic po- 
tentials have been put forward for dark energy, but these 
will not necessarily, and indeed rarely, give rise to peri- 
odic w(a). As one well-studied example the potential 
for a pseudo-Nambu Goldstone boson (PNGB) field can 
be written as V{(j)) — Vb[l + cos((/)//)], clearly periodic, 
where / is a (axion) symmetry energy scale. However, 
unless the field has already rolled through the minimum, 
the relation w{a) is monotonic and indeed can be well 
described by the usual w{a) — wq + Wail ~ a) (see ijlVI 
for further discussion). 

For a damped field, we can estimate how far the field 
has rolled in the history of the universe by writing 

</> ^(TTHo^ VTT^ i?A/p, (1) 

where Mp is the Planck energy and H is the Hubble 
parameter. If w is close to —1, as observations seem to 
indicate, then the field has only traversed 

- </>/i? < Mp, (2) 

and may well not have had time to see nonmonotonicity 
in the potential. That is, if one walks in darkness only a 
few paces one may not realize there are hills and valleys 
in the terrain. 

Therefore, rather than examining nonmonotonic po- 
tentials and calculating periodic behavior in w{a), and 
then investigating those effects on cosmological observ- 
ables, we start directly with a phenomenological, peri- 
odic equation of state (see for some particle physics 



motivation for such an equation of state and 3] for its 
application to crossing w = —1). From the calculation of 
effects on observables we will impose constraints on the 
amplitude and period of dark energy oscillations. 

II. LIMITS ON OSCILLATION 

Our four basic tools for constraining dark energy equa- 
tion of state oscillations will be the distance-redshift rela- 
tion out to redshift z = 2 (as can be accurately measured 
by Type la supernovae), the distance to the cosmic mi- 
crowave background last scattering surface at z = 1089, 
the age of the universe today, and the linear growth factor 
of mass perturbations. These will constrain the models 
sufficiently that we do not need to extrapolate to and 
employ early time constraints from the primordial nucle- 
osynthesis or recombination epochs. However, such early 
universe constraints are important for oscillatory mod- 
els that are not periodic, such as the tracking oscillating 
model of Q. 

Note that despite the double integration leading from 
the equation of state to the distance-redshift relation, 
sufficiently precise distance measurements can provide 
strong constraints on periodicity in 'w(a). Indeed we find 
they are the most incisive probe. 

Since we take a purely phenomenological model for 
if (a), we do not have information on perturbations in 
the dark energy. However, canonical scalar fields gener- 
ically do not have significant perturbations on smaller 
than the Hubble scale, so we will neglect inhomogcneities 
in the dark energy (but see jlllll . Note that coupling of 
dark energy to matter 0, 0, iZ,, can produce both in- 
homogcneities and nonmonotonic equations of state. We 
leave that scenario for future investigation, but gener- 
ally it should cause significant deviations in the growth 
behavior. 

The natural period of the cosmic expansion is given by 
H^^ — {dlna/dt)^^ , so we examine periodicity in units 
of the e-folding scale In a. That is, our ansatz is 

w{a) — wq — Asm{Bhia), (3) 

where wq is the value today. We expect that as the cen- 
tral value, {w) « Wq, becomes less negative, the acceler- 
ation decreases and it will be harder to match the dis- 
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tance and growth relation of, say, a cosmological constant 
w = —1 universe. For example, the undulant universe 
model 13, with an effective (w) — 0, can quickly be ruled 
out on the basis of possessing only 3% of the growth by 
today of the cosmological constant model 7] . 

As the amplitude of oscillations. A, increases we like- 
wise expect clear distinction from a model with constant 
w. As the frequency of oscillations B increases, how- 
ever, the observations, which depend on integrals of the 
equation of state, could have difficulty distinguishing the 
oscillatory model from a model with constant equation 
of state {w). 

To generally describe a sinusoidal equation of state we 
would need four parameters: 

w{a) = Wc - Asm{Bliia + 6), (4) 

where 9 gives the phase of the oscillation today {a — 1) 
and Wc gives the center of the range over which w os- 
cillates (and generally the average value). In Eq. 
and for the rest of this paper, we assume that 9 = so 
that Wc = Wo. This should make it more difficult to dis- 
tinguish an oscillating model from one with a constant 
equation of state w = wq, and so represents a conserva- 
tive approach. (Note that the slinky equation of state 2\ 
has 9 = tt/2, Wc = 0, A = 1, while the undulant model 
further specializes to i? = 1.) 

We explore the three remaining parameters one by one 
to isolate the differing aspects of the oscillations. Of 
course when the amplitude A = then the frequency B 
is not defined as we have reduced to a constant equation 
of state. To actually see oscillatory behavior, we need 
|Blnamin| > 27r, where z^ax = a^ln ~ 1 gives the up- 
per redshift bound on the precision observations, but we 
will consider lower frequencies. For smaller B we might 
simply see a monotonic trend in w over the region of 
observations. Figure plots the redshift of the first full 
period of oscillation. Only for B > 5.72 does a full period 
occur in the range z = — 2. 

For a given wq, we can find the maximum amplitude 
of the oscillation consistent with interpreting the obser- 
vations in terms of a constant, averaged equation of state 
(w) = wq. Thus, larger variations could be distinguished 
as a varying equation of state. As observational con- 
straints we employ future data equivalent to 1% dis- 
tance measures over the redshift range z — — 2 (e.g. 
from SNAP-quality Type la supernovae data 0), 0.7% 
distance measure to the CMB last scattering surface at 
z = 1089 (e.g. from Planck Surveyor CMB temperature 
anisotropy data measure of the current age of 

the universe HqIq, or 2% measure of the linear growth 
factor over z = — 2 (e.g. from SNAP-quality weak 
lensing data). Note that for these probes the cosmologi- 
cal constant {w = —1) case lies respectively 2.3%, 1.0%, 
1.8%, 3.1% away from the w = —0.9 model (keeping the 
matter density fim = 0-28). We find that supernovae 
distances provide the most stringent constraints on the 
oscillation amplitude A and frequency B. 

These results are shown in Fig. [3 for wq = —0.9 and 
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FIG. 1: The frequency B of the oscillatory equation of state 
determines a minimum redshift for seeing a full (or half or 
quarter) period of the oscillations (and hence recognizing a 
distinction from a monotonically varying equation of state). 
The magenta, dotted line denotes z — 2, so for example a full 
oscillation period is only obtained from observations within 
the epoch « = - 2 for B > 5.72. 



three values of frequency B, corresponding to one full 
period (i.e. a true oscillation within the range z = — 2), 
one half period, and B = 1. Larger amplitudes than 
those shown would violate the observational constraints. 
The constraints on A are rather symmetric about zero 
for the small amplitudes allowed; that is, the equation 
of state could equally well decrease into the past by the 
same amount. 

As we take wq less negative, the allowed amplitudes for 
a given frequency decrease. Basically, since the dark en- 
ergy density dies off less quickly for less negative Wq, the 
oscillations have a substantial impact over a larger scale 
factor range, and so constraints are tighter. Upon chang- 
ing Wq we can also adjust the oscillation amplitude A so 
as to keep the same minimum (or maximum) equation of 
state. For example, for a canonical scalar field, w > —1, 
so if we wish to avoid oscillations violating this bound we 
could require wq — A > —I. Figure |31 shows the relative 
constraints on amplitude and frequency when the mini- 
mum is set to the cosmological constant value w — —1. 
For low amplitudes, the oscillations are not discernible 
in the observations so any value of frequency is allowed. 
As the amplitude increases, the frequency must decrease 
in order for the oscillatory equation of state to remain 
looking close to a constant equation of state in the ob- 
servational constraints. 
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FIG. 2: These equation of state curves possess the maximal 
amplitude A allowed for specific frequencies B, given next 
generation observations constraining the magnitude-redshift 
relation to Am < 0.02 relative to the constant equation of 
state model w = wq — —0.9, for z = — 2. Oscillatory 
equations of state with a given frequency must lie closer to 
the dotted, constant equation of state line than the pictured 
curves to match the constraints. Upon flipping the sign of A 
(initially decreasing the equation of state into the past) these 
curves are basically reflected about the w = —0.9 line. 



III. OSCILLATION AT EARLY TIMES 

Oscillating equations of state can act as models con- 
tributing a substantial fraction of the total energy den- 
sity at early times (z ^ 1). These are not really early 
quintessence models, though, in that the equation of state 
is not necessarily providing acceleration, i.e. w < —1/3. 
The energy density of the oscillating model is 

Posc(a)=Posc(l)a-^(^+'""^e(3^/^)[i— (5) 

For the limiting cases of Figs. I2I3I the contribution to 
the total energy density at z = 2 is within 20% of what a 
plain w = —0.9 model would give, and the energy density 
fraction at z = 1089 remains at the 10~* level. The 
slinky model manages to have high early energy density - 
possibly even greater than the matter density - by taking 
the central value Wc = 0. Thus, early energy density is 
achieved at the price of not maintaining dark energy, i.e. 
acceleration. 

Another caution is that early presence of substantial 
energy density not in the matter component can obvi- 
ate the usual matter dominated epoch. While the con- 
sequences of this for structure formation are recognized 



FIG. 3; For a frequency B the curve defines the maximal 
amplitude allowed, given next generation observations con- 
straining the magnitude-redshift relation to Am < 0.02 for 
z = — 2 relative to the constant equation of state model 
w — Wq = —1 + A. Such a model has an equation of state 
reaching a minimum value of —1, the cosmological constant 
value. 



in the suppression of the linear growth factor, there is a 
more subtle effect. The second order differential equa- 
tion for the growth factor g = {5pm/Pm)/a requires two 
boundary conditions, canoni cally given by g(a ^ 1) = 1 
and dg/da{a ^ 1) = (cf. [IJ). However, these arise 
from the usual density growth behavior Spm/ Pm ^ a for 
a matter dominated universe. In numerical computation 
of the growth equation in an oscillating model, we find 
a sensitive interplay between the initial scale factor 
and the fraction of energy density at that scale factor 
in components other than matter, such that the bound- 
ary condition is unstable. Consequently, for slinky mod- 
els, numerical instability of the growth factor arises for 
B > 0.2. This may be apparent in the nonmonotonicity 
of the growth shown in Fig. 11 of 9]. Proper solution 
in this regime requires solution of the full set of coupled 
matter, scalar field, and radiation perturbation equations 
so as to give the proper initial conditions (velocity) to the 
matter growth. Note that some initial velocity dg/da > 
can actually help offset the suppression of the growth fac- 
tor. 

Models, e.g. where Wc ~ —1, that do have a true mat- 
ter dominated epoch, do not suffer from these complica- 
tions. Interestingly, if Wc = wq ^ —1, growth data alone 
allows two distinct regions in parameter space. For low 
frequencies B, the equation of state is nearly constant 
and so the growth mimics that of the constant w model. 
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(If Wc = —1, the oscillating model in the limit S — > 
approaches the cosmological constant.) As B increases, 
the growth is altered due to the change in the Hubble 
parameter relative to the constant w model. The tension 
between growth and solution of the coincidence problem 
(that acceleration is observed near today) is nicely high- 
lighted by p. For very high frequencies B, the growth 
once again approaches the constant w case (cf. Eq. [Sjl . 
This holds even if the amplitude of oscillations A ^ 1. 

A high oscillation frequency, though, may be expected 
to be accompanied by spatial inhomogeneities in the 
dark energy due to the (relatively) large effective mass 
of the field. This will then alter the matter perturba- 
tion growth. Regardless, use of precision distance data 
removes the high frequency fit region, for large ampli- 
tudes. For small amplitudes any frequency B is allowed 
by the distance data, as shown in Fig. |2 

IV. OTHER OSCILLATIONS 

We briefly return to the issue of oscillations of the po- 
tential, or the field motion in the potential, rather than 
oscillations of the equation of state directly. As men- 
tioned in ^ the PNGB model has a periodic potential. 
For accelerating behavior, however, the field must not 
have undergone many oscillations, else the equation of 
state would be near zero like the axion case. 

Caldwell & Linder |l3j| examined the dynamics and 
found the behavior in the w' — w phase plane to be nearly 
linear: w' ~ F(l + w), where w' = dw/dlna. The coef- 
ficient F is proportional to the inverse of the symmetry 
scale /. They showed this was a thawing model, where 
the field starts frozen and hence looking like a cosmolog- 
ical constant w = —1, before the field begins to evolve 
as it comes to dominate the cosmic energy density. The 
solutions for the equation of state and energy density for 
such dynamics are 

w{a) = -l + (l + wo)a^ (6) 
P«,(a) = p.(l)e(3/^)(i+-»)(i-'^") , (7) 

which are manifestly non-oscillatory. As mentioned in 
^ this looks like the standard varying equation of state 
'w{a) ^ wq + Wa{l — a) when _F « 1, which holds for 
natural values / ^ 0{Mp). 

Another possibility is for the potential not to be peri- 
odic but the field motion is. Consider oscillations around 
a minimum of the potential. Turner [31 showed that for 
oscillation periods much smaller than the Hubble time 



(large frequency B in our notation), the averaged equa- 
tion of state becomes 



for potentials y ~ 0" with n even. For potentials looking 
like quadratic (quartic) potentials near their minimum, 
this gives (w) = (1/3) - behavior like matter (radi- 
ation). To provide acceleration, one would need n <C 1 
(conversely, very steep potentials give behavior like shear 
energy, (w) =1). If the frequency of oscillations is com- 
parable or larger than the Hubble time (so the energy 
density changes significantly during the oscillations), the 
equation of state does not obey Eq. but in this case 
we would not discern oscillations in the data. 

Thus, it is not easy to obtain oscillations in the equa- 
tion of state without fine tuning the potential to have 
particular nonmonotonic periodicities, as in the slinky 
model. We can analyze the phase plane dynamics of our 
oscillatory model, Eq. quite simply (also see Fig. 9 
of 9]). It satisfies the equation 

SO the dynamics is given by an ellipse in the phase 
space with center at w = wq, w' = 0, semimajor axis 
A, semiminor axis AB, and eccentricity y/l — B^ (or if 
B > 1, switch major and minor, and the eccentricity is 
Vl — B~'^). Note that B = 1, where the characteristic 
oscillation scale equals the Hubble e-fold scale, gives a 
circle. 



V. CONCLUSION 

Precision distance-redshift measurements over the 
range z — 0—2 possess power in discriminating oscillating 
equations of state from constant ones, despite the double 
integral relation. Other probes are not as constraining, 
but growth measurements can provide important limits 
on the early time behavior if the dark energy affects the 
degree of matter domination. Growth must be calcu- 
lated carefully, however, since high amplitudes affect the 
growth boundary conditions and numerical stability, and 
high frequencies can lead to spatial inhomogeneities in 
the dark energy. While oscillating models offer one idea 
for solving the coincidence problem (since acceleration 
occurs periodically), next generation observations will be 
able to strongly constrain such possibilities. 
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